Analytic results on long-distance entanglement mediated by gapped spin chains 
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Recent numerical results showed that spin chains are able to produce long-distance entanglement 
(LDE). We develop a formalism that allows the computation of LDE for weakly interacting probes 
with gapped many-body systems. At zero temperature, a dc response function determines the ability 
of the physical system to generate genuine quantum correlations between the probes. We show that 
the biquadratic Heisenberg spin-1 chain is able to produce LDE in the thermodynamic limit and 
that the finite antiferromagnetic Heisenberg chain maximally entangles two spin- 1/2 probes very far 
apart. These results support the current perspective of using quantum spin chains as entanglers or 
quantum channels in quantum information devices. 

PACS numbers: 



The quantum information (QI) approach of regarding 
entanglement as a information resource pj stimulated im- 
portant developments in its characterization and in ways 
to measure it. On the other hand, in condensed mat- 
ter physics, quantum correlations have long been recog- 
nized as an essential ingredient in many low-temperature 
phases. These facts lead to a considerable amount of 
work on the characterization of the entanglement prop- 
erties of many-body systems at zero temperatures (for 
reviews see [2]), particularly near quantumphase transi- 
tions and also at finite temperatures 0, 0, 0] . 

Feasible mechanisms of entanglement extraction from 
real solid state [6| and their ability to transfer entangle- 
ment between distant parties [3| are of crucial importance 
for the implementation of QI protocols, such as telepor- 
tation or superdense coding. In systems with short-range 
interactions, though, entanglement between two particles 
usually decays quickly with distance between them 
However, Campos Venuti et al |3], have recently found, in 
numerical density matrix renormalization (DMRG) stud- 
ies, that certain correlated spin chains are able to estab- 
lish long-distance entanglement (LDE) between probes 
to which they couple, without the need of an optimal 
measurement strategy onto the rest of the spins Q • This 
naturally raises the question of which classes of strongly 
correlated systems are able to produce LDE. 

In this paper, we present a quantitative description of 
the effective Hamiltonian of interaction between probes 
weakly coupled to gapped many-body systems. This al- 
lows us to confirm some of the numerical evidence of LDE 
found by DMRG @ and to obtain new results concerning 
two important spin chain models. Although our formal- 
ism can be applied to general gapped many-body sys- 
tems, here we focus on one-dimensional spin chains. In 
the spirit of @, two probes, a and b, interact with the 
spin chain, locally, through sites m and n, respectively. 
The Hamiltonian of the system reads H = Hq + V mi „, 
where H is the full many-body Hamiltonian of the spin 
chain and V m ,n = V a ,m+Vh,n describes the interaction be- 
tween the probes and the spin chain. We will show that, 
as long as the probes interact weakly with the many- 
body system, the ground state (g.s.) of the full system 
may display LDE between the probes, i.e., E(p a b) > 



when \m— n\ is of the order of the system size. The oppo- 
site limit, i.e., strong interactions, will cause the probes 
to develop robust correlations with the site they inter- 
act inhibiting entanglement with them @|. This arises 
from a constraint on the correlations between different 
sub-systems known as the monogamy of entanglement 
[H]. In order to maximize the spin chain potential to en- 
tangle the probes, we require that J p <C J, where J p 
is the interaction strength between the probes (a, b) and 
the spins at m and n, and J is a typical energy scale 
for the spin system (for instance, a nearest-neighbor ex- 
change interaction) . For J p = the state of the probes 
becomes totally uncorrelated and the g.s. of the entire 
system becomes (d x <i)-fold-degenerate. In this case we 
may write \ip) = \ip ) ® \x a ) <8> \Xb), where \ip ) is the g.s. 
of the spin chain (assumed non-degenerate) and |x 7 ) the 
state of the probe 7. The role of the interaction ( J p > 0) 
is to lift this degeneracy causing the probes to develop 
correlations. For weak coupling J p , the effective Hamil- 
tonian in this d x d low-energy subspace, obtained as 
discussed below, will determine the dynamics and corre- 
lations of the probes. For the special case of spin one-half 
probes (d = 2), the negativity [13], or any other equiva- 
lent entanglement monotone, can be used to quantify the 
entanglement. 

The effective Hamiltonian. In a very general way we 
can write the local interaction between the probes and 
the corresponding sites in the many-body system in the 
following manner: 

V = j2l a a O« n ®A a ®l b + 1 b p«®l a ®B a , (1) 

Q = l 

where A a (B a ) denotes an operator acting on the Hilbert 
space of the probe a(b) and l a (b) the corresponding iden- 
tity operators. The many-body system operators on 
sites to are represented by and 7« is the cou- 
pling strength. The projector onto the states with unper- 
turbed energy £0 = (ip\Ho\i/>) isVo = \ifo) (V>o| ® la ® lb- 
Let Vk (k > 0) be the projector onto the subspace 
of energy E k > E a , so that 1 = Vq + 52 k>0 Vk, 
Using the standard canonical transformation formalism 
[m one can determine probes g.s. by diagonalizing 
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an effective Hamiltonian in the subspace spanned by 
Vq, namely, H eff = V (H + V)V + |Vo)<V>o| ® # (afc) 
[121 ] . This a familiar concept that finds many appli- 
cations in condensed matter physics, such as, for in- 
stance, in the derivation of the Ruderman-Kittel-Kasuya- 
Yosida magnetic interaction between local moments in a 
metal [13]. The coupling between the probes is given by 
= -J2k>o( UP kU)o(E k - E )-\ where the aver- 
age is taken with respect to the g.s. of the spin chain 
(UP k U) = (^ \UP k U\%) and U := V - (V) . Note 
that, by definition, (U)o = VqUVo = 0. Entangle- 
ment between the probes arises from ij( ab ) since it con- 
tains nonlocal terms such as U a , m 'PkUb y n [HI- The probe 
Hamiltonian can be transformed by straightforward ma- 
nipulations into an explicit form involving time depen- 
dent correlation functions of the spin chain. A similar 
procedure is used to express cross sections of scattering 
by many-body systems in terms of its correlation func- 
tions [23J. We obtain (we set H = 1, see [15( for the 
derivation) 

H(ab) = ~hj lt\ dt ( U (W)oe iEt . 

We now introduce the explicit form of U to arrive at 
the desired result = + + H$> + H^ b) , 

(see [l(| for a comment on the local terms). Defin- 
ing the two-body connected correlation {0" n (t)0^) c — 
(^ \O^(t)O^ ) - (^o|0£(<)|^oWo|Qg|lfo) the term 
coupling the two probes yields 

p 

h nl = Yl l a a ll{C ma , n p + C nP , ma )A a ® B0(2) 

a, 13=1 

C„ & = - J dte~° l t l S i 5 n(t)<O«(t)O^(0)) c .(3) 

The coupling between the probes can be expressed 
in terms of the response function Xma;np(t) = 
-i([0° (*), O£])0(t), where 9{t) is the Heaviside step 
function. Using the Lehman representation at T = 
one can show that 

where Xma-nf3(to) is the time Fourier transform of 
Xma-np{t)- The effective Hamiltonian Eq. (f2]) lifts the 
degeneracy of the GS level of the uncoupled system 
(J p = 0). As long as the couplings appearing in -ffjvx* are 
small compared to typical energy scales of the spin chain, 
such as the gap to first excited state, A, the low-energy 
physics of this system, E <C A, with no real excitations of 
the spin chain, will be well described by . This con- 
dition limits the strength of the chain-probe interaction, 
but is shown by numerical results to be the appropriate 
limit to maximize LDE p|. 

In the remainder of the paper we will compute the LDE 
for two rotational invariant spin chains: the finite Heisen- 
berg spin-1/2 chain in zero field using the exact results 



from bosonization theory; a specific spin-1 Heisenberg 
chain with biquadratic interactions by means of a ap- 
proximation scheme for its spectrum. It is useful to write 
Eq. (Q]) in terms of spin operators S m for the spin chain 
and r (M for the probes. Considering that the probes cou- 
ple with the spin chain via an Heisenberg interaction, the 
most common situation, V m , n — JaS m ■ f a + J] } S n ■ t%, the 
connection with the previous notation becomes straight- 
forward: 0° = S°, A a (B a ) = T« (b) and 7 a (fc) = J a(b) . 
Eq. ([2]) becomes simply, 

Heff = JabTa • U, (4) 

where J a b = J a JbXma;na{0) and a = x, y or z. 

The finite antiferromagnetic Heisenberg spin-1/2 
chain. The isotropic antiferromagnetic Heisenberg model 
reads 

L-l 

H := Si ■ S i+1 . (5) 

Our formalism only applies to the finite chain which has 
gapped excitations. To calculate its time-dependent cor- 
relation functions we will use the conformal invariance 
of the critical infinite chain (L — > oo). We show below 
that its time-dependent correlations are enough to ex- 
tract the effective coupling J a b for the finite chain. It is 
clear that the effective Hamiltonian Eq. $2$ will preserve 
the full SU(2) symmetry of the interaction Hamiltonian 
H + V mt n, i.e. no local terms will give contribution to 
H e f f . Assuming that the probes couple to the spin chain 
with the same strength (J a = Jb = J p ), H e ff takes the 
very compact form 

Heff = JpXma;na(0)f a ■ f b . (6) 

Hence, whenever Xma;na(Q) > the g.s. of the probes 
is a singlet-state displaying maximal LDE - Our 
computation of Xma-na(0) will rely on general bosoniza- 
tion results for correlations of spin-1/2 chains (see [l8| 
for a review) and the conformal invariance of the crit- 
ical chain. The dominant long-distance correlations of 
the g.s. of Hamiltonian Eq. |(5]) oscillate with a ir 
phase change between neighbor spins. It is therefore 
useful to define the retarded Green function for the 
staggered magnetization Mf := (— l) J '5f, G^ n (t) := 
i([M%,M%(t)])d(t). The response function in Eq. © 
is Xmocynccit) = (-1) |m_ "' C (t) . The retarded Green 
function is obtained from the corresponding Matsubara 
Green function, G(x,t) := {f T M%{x m ,T)M%(x n ,Q)), 
with imaginary time r S [— (3, (3} and where, in the limit 
\%m — x n \ ^> 1, we may replace x m — x n by a contin- 
uum variable x. The Matsubara Green function for the 
infinite chain reads 

G(x, t) = A\vpr + ix\~ ; 

A is an amplitude and vp the Fermi velocity of the spinon 
excitations. This result implies that the infinite chain 
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Figure 1: The conformal transformation w(z) = u + ir maps 
every point (vfT, x) in the plane into the strip geometry (u £ 
] — oo, oo[, r 6 [-L/2, L/2]) with periodic boundary conditions 
along the v direction. 



has a divergent Xmmnai®)', this is a direct consequence 
of a zero gap and a signal of the critical nature of the 
spin chain at T = 0. Nevertheless, the finite chain is 
gapped, has a finite Xma-na(0) , and the above result can 
be used to calculate its correlation functions, precisely 
because the critical nature of the infinite chain implies 
that iV-point correlation functions in different geometries 
are related by conformal transformations [20]. The map- 
ping of the infinite chain to the finite chain is achieved 
by the following analytic transformation (see Fig. 1), 
w = In z = u + ir 1 where z — vft + ix. 

Using the transformation law for conformal invariant 
theories [13] the Matsubara Green function for the fi- 
nite antiferromagnetic Heisenberg chain with periodic 
boundary conditions in the spatial coordinate r reads (2l| 
G cy i(r,u) = 27rf (2cosh(27rw/£) - 2 cos(27rr/L))"-?.The 
analytic continuation to real time is made by Wick ro- 
tation u — > ivpt + + sgn(t) and the corresponding re- 
tarded Green function defined in the cylinder G^ yl (r,t) 
can be computed from the time-ordered Green function 
G R (r 7 t) = -20(*)Q [G(r,t)](see [H]). Setting the branch 
cut of the logarithm in the negative real axis we find 



G*(r,t) = 2n 



A 6(t)6(F(r, t))sign(sin(2wt/L)) 



where F(r,t) = 2cos(27rr/L) - 2cos(2vrf/£). The re- 
sponse function at zero frequency is then given by 
Xma-MO) = (-1)1™-™! f °° dtG R {r,t)exp{-0+t). Here 
we only state the result, 



X 



ma.;na 



(o) = c„ 



dy 



y/ir - 1 



-\/cos(27rr/L) — cos(y) 



with C mn = (-l) |m ~™U/(2u f ). Figure 2 shows the plot 
of the absolute value of the response function at zero 
frequency confirming the existence of LDE for a wide 
range of values of r/L. Note that x mtt;ntt (0) diverges 
logarithmically at the origin. Our perturbative approach 
cannot be applied unless JpXma-na{0) € A~ J/L, and 
will fail in the thermodynamic limit (L — > oo) for fixed 
r. The numerical results of Ref. Q show probes almost 
completely entangled only for small values of J p ~ 0.1, 
for a finite chain L = 26. This value is well estimated 
by the limit of validity of our perturbative approach, for 
r/L ~ O(l), namely, J p J/^/L. These results, in the 
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Figure 2: The absolute value of the response function at zero 
frequency for the finite antiferromagnetic Heisenberg chain. 
We have assumed r > 1 so that the results from bosonization 
theory are accurate. 



light of our analysis, strongly suggest that the conditions 
for LDE are coincident with the conditions for validity of 
the perturbative approach. Weakly coupled probes get 
maximally entangled by the effective antiferromagnetic 
interaction mediated by the spin chain. 

The AKLT model. The Heisenberg spin one-chain with 
biquadratic interactions reads 



H := 



N-l 

E 



Si ■ Si+i + [3{Si ■ Si+i)' 



This model admits an exact solution for (3 = 1/3 which 
is known as the Affleck-Kennedy-Lieb-Tasaki (AKLT) 
point [22]. A picture of the g.s. is given by the so-called 
valence-bond-solid (VBS). Each spin-1 is represented by 
a couple of spins one-half, as long as the antisymmetric 
state is projected out. The VBS state is constructed by 
forming short-ranged singlets between nearest spin- 1/2 
and then symmetrizing local pairs to get back S = 1 
states. In the thermodynamic limit, the static correla- 
tions are very short-ranged [£,aklt — l/ln(3) = 0.9] [22T |. 
For this reason we may ask whether two probes are able 
to get entangled by interaction mediated by the spin-1 
chain. We cannot make an exact computation of LDE as 
in the Heisenberg model, since the exact dynamical cor- 
relations are not known even for l arg e distances. How- 
ever, as suggested by Arovas et al [2J|, we can apply the 
single-mode approximation (SMA)used to deduce the 
phonon-roton curve in liquid 4He [25|, in order to study 
the excitations in this model. This is done by assuming 
that a excited state at wave vector q is given by 

\q) = S^o) =N- 1 / 2 J2 etqr ' S ^o), 

i 

where IV'o) is the exact g.s. of the AKLT model. 
Within the SMA the dynamical structure factor S a/3 
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is related with the static structure factor defined as 
s a < 3 {q) = {ip \S-qS^\fo) in the simple way S{q,w) S 
s(q)S(w — w q ). In [241] it was shown that, w q = E q — E = 
5(5 + 3cosg)/27 and that s(q) = (10/27) (1 - cosq)/w q . 
The knowledge of the dynamical structure factor allows 
us to compute the effective couplings of Eq. by inverse 
Fourier transform. For the AKLT model we obtain 

f°° dt r b 

Xm+m;ma(0) = / I dqcos(qr) sm(w q t)(a -\ ) 

JO 27T J-K W q 

with a = —2/3, b = 80/81. These integrals may be 
done by defining f(t) = sm(w q t)8(t) and noting that 
7 = f^dtsmfaq^O^) = f(w = 0) where f(w) is 

the Fourier transform of f(t). We obtain 2f(w) = 
(w + w q — i0 + )~ 1 — (w — w q — i0 + )~ 1 and 7 = l/w q . 
The remaining integral is done by extending the inte- 
grand to the complex plain and computing the residues. 
This yields, J ab = -§J*(-l) r (l + §r)e~*I#Z7. The 
sign of the interaction mediated by the AKLT spin chain 
changes according to the distance between the probes. 
This comes from the fact that the static correlations in 
this spin chain have a similar alternation. Therefore at 
T = the probes get entangled whenever their distance 
corresponds to a odd number of sites. 

The effect of finite temperature and final comments. If 



the temperature is such that fceT <C A, we do not expect 
real excitations of the spin chain to be present. Only the 
subspace of states described by H e f / will be populated 
and we may calculate the correlations between the probes 
using p ab = e-P&ofS /Tr[e-P A °ff] with 0- 1 = k B T. This 
defines a temperature, 1/0*, above which entanglement 
disappears. For an antiferromagnetic H e fj the compu- 
tation of the negativity yields 0*J mn ~ 0.27. Loosely 
speaking, the probes will be entangled whenever the tem- 
perature is smaller than the effective coupling between 
the probes. This formalism has also been used to dis- 
cuss qubit teleportation and state transfer across spin 
chains j§]. In conclusion, we have expressed the capac- 
ity of a gapped many-body system as an entangler of 
weakly coupled probes at arbitrary distances in terms 
of a zero temperature response function. We exempli- 
fied this formalism by calculating this function for two 
quantum spin chains, shedding light on recent numeri- 
cal results on LDE. In this context, our results strongly 
suggest that the main mechanism of LDE mediated by 
gapped spin chains is the existence of dominant antiferro- 
magnetic correlations. Acknowledgments. We gratefully 
acknowledge very enlightening discussions with J. Pene- 
dones. A.F. was supported by FCT (Portugal) through 
PRAXIS Grant No. SFRH/BD/18292/04. 
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